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I. INTRODUCTION 



A microscopic multicluster model is successful for understanding the structure of light 
nuclei [lj . One example recently applied to halo nuclei is the description of lithium isotopes 
7 ~ n Li in terms of an a + t + xn model including (2 + x) clusters, where x denotes 
the number of neutrons, i.e., x = 0, 1, ... ,4 for 7_11 Li. With the use of an effective two- 
nucleon potential, this approach reproduces consistently the ground state energies of these 
nuclei together with their proton and neutron radii. In particular, a recent experiment [3] 
has confirmed that the isotope dependence of the charge radii is best reproduced by the 
stochastic variational calculation |4J on the basis of this multicluster model. 

In the microscopic cluster model, the structures of clusters are described in terms of 
appropriate wave functions, usually in the simplest configurations of the shell model, and 
an antisymmetry requirement on the total wave function is taken into account exactly. The 
multicluster model is thus a microscopic many-body theory based on an effective nucleon- 
nucleon potential. One must, however, compute fairly complicated matrix elements for those 
non-orthogonal basis functions used in the multicluster model, which makes it hard to apply 
to larger nuclei. This motivates another treatment of the multicluster system, that is, a 
macroscopic cluster model. 

In a macroscopic cluster model, the clusters are treated as structureless point particles. In 
this case, cluster-cluster interactions are usually represented by appropriate local potentials 
which are, more or less, phenomenologically determined so as to reproduce the relevant data 
on the cluster-cluster systems. One of the simplest and well-investigated examples is a 3a 
model for 12 C. 

Many calculations based on the microscopic 3a model show that it reproduces the general 
features of the low- lying spectrum and the spectroscopic properties of 12 C p, |6(. Thus we 
would expect that the macroscopic version of the 3a system could give results of similar 
quality if it would be a realistic substitute of the microscopic model. Since the a particle 
is a tightly bound system, it is natural to conceive that this expectation may come true. 
Two kinds of local a — a potential are often employed in the macroscopic 3a model. One 
is an angular momentum dependent potential such as Ali-Bodmer (AB) potential Q, or 
that proposed in Ref . [|[ . These potentials are repulsive at the short a — a distance in 
the s and d waves. The other is a deep, angular momentum independent potential such 
as Buck-Friedrich- Wheatley (BFW) potential j^, • Both potentials reproduce the phase 
shifts of low-energy a — a elastic scatterings well, but they differ in the prediction about 
a 2a bound state. The AB potential produces no bound state, consistently with the fact 
that the 2a system has no bound states. Contrary to this, the BFW potential is so deep 
that it produces two bound states (0s and Is) in the s wave and one bound state (OeQ in 
the d wave. These bound states are considered redundant or forbidden states from which a 
physically acceptable a— a relative motion function must be free. The removal of unphysical 
bound states has been done using either the projection technique with an orthogonalizing 
pseudopotential Ar 11 1 or a supersymmetric transform of the original potential 12 . It was 



shown in Ref. [12j that the supersymmetric partner potential of the BFW potential is similar 
to the AB potential. 

How well does the macroscopic model with the AB or BFW potential reproduce the 
results of the underlying microscopic calculation for 12 C? In this respect, the macroscopic 
calculation gives quite poor results Til 3 IE 13, e -S-> the ground state energy appears just 



below the 3a threshold, which is far from the experimental value, —7.27 MeV. Moreover, the 
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calculated ground state hardly has the properties suitable for the 12 C ground state. It is also 
confirmed that the projection technique gives highly unstable results which depend on the 
projection constant A and on how accurately the forbidden states are solved as well [l3L Eif. 
A convergent result is obtained only by letting A sufficiently large. These results are quite 
embarrassing considering the success of the microscopic 3a model. 

The purpose of the present study is to clarify the origin of the enigmatic behavior of the 
projection technique as well as the reason why its ground state obtained in the macroscopic 
3a model is far from the physical ground state. To answer these questions, we carry out a 
direct method of removing the forbidden states by diagonalizing the projection operator T. 
We will show, using an example of the BF W potential, that the enigma can be understood 
by the existence of particular 3a states, which are excluded from the 3a function space in 
the projection technique, but, in our opinion, should be included as allowed states. This is 
demonstrated by defining both of allowed and forbidden states in translationally invariant 
harmonic-oscillator basis functions. 

In the microscopic 3a model, the a particle is usually described in terms of a harmonic- 
oscillator (Os) 4 configuration. Then it is natural to define allowed or forbidden states re- 
ferring to the Pauli exclusion principle. Such allowed states are obtained explicitly for 
the 3a system through the norm kernel of the microscopic 3a wave function [17]]. A sim- 
pler approximation is to remove the 2a Pauli-forbidden states, namely the Os, Is and Od 
harmonic-oscillator states, from any pair of the a — a relative motion 0, 0] • We call the 
resulting states "Pauli-allowed" states, and distinguish them from those defined according 
to the redundant bound states in the BFW potential. 

In sect. HU we discuss the methods of removing the forbidden states in the 3a system. We 
stress the utility of the harmonic-oscillator basis functions, which makes it straightforward 
to construct the allowed states. In sect. IHIl we discuss the Pauli-allowed states which have 
direct relationship to the microscopic 3a model, where the a particle is described in terms of 
a harmonic-oscillator (Os) 4 configuration. In sect. IIVI we show how the enigmatic behavior 
of the projection technique can be understood. In sect. W\ we discuss such an a-a potential 
that is suited for the macroscopic calculations. A summary is given in sect. IVII Some of the 
lowest Pauli-allowed states are given in Appendix. 

II. SOLUTION WITH THE REMOVAL OF FORBIDDEN STATES 

The Hamiltonian for the 3a system is in the form 



where m a is the mass of the a particle and the kinetic energy of the center of mass motion 
is subtracted. Vi is the a— a potential acting between (jk) pair. Here we use the convention 
of (jk) = (23) for i = 1, (jk) = (31) for i = 2, and (jk) = (12) for i = 3. We want to solve 




(1) 



(2) 



with the constraint 



(p/fo) I *) = (i = l,2, 



3) 



(3) 
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where tpf stands for a particular forbidden state. When we have several forbidden states, 
the above condition must be applied to all of them. The a — a relative distance vector is 
given by r*j = Rj — Rk, where Ri is the center of mass coordinate of the zth a particle. 
Another relative coordinate is denoted by p { = (Rj + Rk)/2 — Ri. The wave function \1/ 
must be totally symmetric with respect to the interchange of iZj's. 

As in the BFW potential, we assume that there exist three forbidden states, Os, Is and 
Od. The function ipj in Eq. (jHJ) thus represents these three states. To solve Eq. (J2J) with the 
condition ©, we define an operator 



r = £ r - (4) 

i=i 

where 

r^r^ + r^ + r^od) (5) 

with 

1 

Ti(n£) = 2J I ¥n£m(ri))((Pnem(ri) I . (6) 

m=—£ 

The operator Tj is a projector onto the forbidden states for the relative motion between 
the a particles of (Jk) pair, but we must note that, for i ^ j, it satisfies the relation 

r,r., / r ; r, / 0. 

We consider the eigenvalue problem of F 

rxS^TxS, (7) 

and obtain those eigenf unctions which are totally symmetric, i.e., [3] symmetry. The eigen- 
value 7 is non-negative. This is seen as follows: 

7 = (xS I r I xf M ) 

= 3 (xffi I Ti | X 1 !m) (< = 1,2,3) 

= ^<X^\<Pf(ri))(<Pf(ri)\x^)- (8) 

Here we used the symmetry of x^lm an d noted that the matrix element in the last line is 
non-negative, which follows from the fact that ((pf(ri) \ Xlm) ^ s a function of p iy say, v(p i ), 
and thus (Xlm I V 9 /( r i))(v 9 /( r j) I Xlm) = ( v I v ) can no ^ ^e negative. Clearly 7 = is 
equivalent to v — | Xlm) = (* = 1)2,3), namely, the eigenfunction with 7 = 

contains no forbidden components, so it satisfies the condition (jSJ). We can thus expand the 
solution of Eq. (J2J) in terms of the eigenfunctions belonging to 7 = 0, 

* = E^xS 7 m°, (9) 
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where p is introduced to distinguish the orthonormal eigenfunctions with 7 = 0. The 
coefficients C p and the energy E are determined from the secular equation 



j:(x^M°\H-E\ X ^°)C p/ = 0. (10) 
P > 

The procedure explained above is a standard way to obtain a solution which is free from 
any forbidden components. One, however, has to note that to obtain vanishing eigenvalues 
of Eq. (J7J) demands a calculation of high accuracy 

Another method of solution equivalent to the above direct method is to add the orthog- 
onalizing pseudopotential to the Hamiltonian 

H = H + XT (11) 

and solve the equation for \I/ ( A) of [3] symmetry 

HV{\) = E{X)V{\) (12) 

as a function of the projection constant A. The solution for \I/(A) may be expanded as in 
Eq. © but then all the eigenfunctions with both 7 = and 7 7^ must be included. The 
desired solution is such that E(X) and ^(A) are stable for sufficiently large A. To understand 
this, we note that the energy E(X) reads as follows: 

E(X) = (tf (A) I H I (A)) + A(tf (A) I T I (A)). (13) 

The first term on the right-hand side of Eq. (|13|) is expected to change moderately as a 
function of A, but the second term increases rapidly with increasing A in such a way that 
^(A) has a significant overlap with the forbidden states because in that case the probability 
of finding the forbidden components in ^(A), (^(A) | T | ^(A)), becomes a positive, non- 
vanishing value. In contrast to this, E(X) can be stable for sufficiently large A, provided that 
^(A) has vanishingly small or no overlap with the forbidden states. The advantage of this 
pseudopotential method is that one does not need to solve Eq. (J7J) to construct the allowed 
space spanned by the solutions with 7 = 0, but only needs to diagonalize the Hamiltonian 
with the pseudopotential. To obtain a stable solution, however, one has to provide a function 
space large enough to remove the forbidden states. 

It is important for our purpose to construct a complete set of [3] symmetry for the 3a 
system. For example, the eigenvalue problem (j7J) of the operator T can be solved in such a 
set. The construction of the complete set can be performed using translationally invariant 
harmonic-oscillator functions 

K(nX(Pi)]LM, (14) 

where ip q lm is the oscillator function with the number of oscillator quanta q, the orbital 
angular momentum i and its z component m. (q—£ must be a non-negative even number.) 
The symbol [• • - ]lm indicates the angular momentum coupling. The oscillator parameters 
of the functions ^^(fi) and ^ q i 2m2 {P\) are 2u and |z/, respectively. Here v = (m N u/2h) is 
the oscillator parameter of the single-nucleon wave function in the harmonic oscillator shell 
model, where is the nucleon mass. We choose v = 0.26 fm~ 2 because it is an appropriate 
value to reproduce the charge radius of the a particle using the (0s) 4 wave function [I?}. 
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One may choose other values of u, but then the diagonalization of the operator T and the 
Hamiltonian, Eqs. (J7|) and ()10j). would require larger dimension. Though an elegant method 
of constructing the states with definite permutational symmetry is given in Ref. 20], we 
follow a simple procedure of diagonalizing the symmetrizer 

5 = i(l + (1,2,3) + (1,3,2)), (15) 

where (1,2,3) and (1,3,2) are cyclic permutations. The states with [3] symmetry can be 
obtained as those with unit eigenvalue of S. The symmetrizer S conserves the number of 
total oscillator quanta Q = qi + q 2 , so the diagonalization of S can be done separately 
for each Q. In this diagonalization one needs the matrix element of the permutation P 
((1,2,3) and (1,3,2)) in the basis (JHJ), which can be easily evaluated using a generalized 
Talmi-Moshinsky bracket |2oL l2ll ] . 



III. PAULI- ALLOWED STATES IN HARMONIC-OSCILLATOR BASIS 

So far we have assumed that the forbidden states <ff(r) are the bound states of the a— a 
potential, which is considered natural in several studies [ill EH HE EE HI- A different 



viewpoint is that the forbidden or redundant states should be defined according to the 
underlying microscopic model, that is, they are determined solely from the wave function 
of the a particle, not from the a — a potential employed. For the a particle described with 
the (0s) 4 harmonic-oscillator shell-model with the oscillator constant u, it is well-known 
that three Pauli- forbidden states appear for the a — a relative motion. They have the same 
quantum numbers, 0s, Is and Od, as those of the BFW potential, but their wave functions 
are simply given by the harmonic-oscillator functions ^L,((/ = 0, 2). This approach is in fact 
a multicluster version of the orthogonality condition model 0, 22|. An a— a potential used 



in this case is often taken to be close to a double folding potential [2c 

It is useful to define an operator r HO , an analogue to the operator ([%[). in the harmonic- 
oscillator basis, where Ti{nt) of Eq. © is replaced by 

t 

rf°K) = £ I V-LMWLfo) I (g = 2n + £). (16) 

m=— I 

The eigenvalue problem for r HO , the analogue of Eq. (J7J), can be solved in the basis (|14|). 
The eigenf unctions can be classified by Q, an SU(3) irreducible representation (A/x) and an 
additional quantum number k [3, [3] (k serves to distinguish a multiple occurrence of the 
same (A/x) states, and it is otherwise suppressed). That is, 

xfS (Xfl)K \ (17) 

and the eigenvalues 7 are either zero (Pauli-allowed) or non-zero values of order unity (Pauli 
forbidden). The SU(3) coupled basis is easily constructed by use of SU(3) coefficients 



[^(r)^(p))^ = £((910)4, (g 2 0)i 2 || (A/i)L) [V£(r)Vg(p)]LM. (18) 
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Here (A/i) can in general take the representations of (qi + q 2 , 0), (q\ + q 2 — 2, 1), (gi + q 2 — 
4, 2), ... , (q> — q<, q<), where q > = max(gi, q 2 ) and g< = min(gi, q 2 ), but is limited to those 
which contain the angular momentum L. For L = 0, both A and // must be even, and a 
concise expression for the SU(3) coefficients with L = is given in Ref. [iif . 

For Q < 8, all the [3]-symmetry basis states are Pauli-forbidden. For Q = 8 we have only 
one Pauli-allowed state with (A//) = (04). For Q = 10, two Pauli-allowed states, (24) and 
(62), appear. We give in Appendix the Pauli-allowed states -^-J^^^ 7-0 f or Q = 8 — 14 in 
the SU(3) coupled basis of Eq. (fTS|) . 



IV. CASE OF THE BFW POTENTIAL 



The BFW potential reproduces very well the a — a phase shifts of the s, d and g waves 
for -Ec.m. < 17 MeV 0. It is a local, angular momentum independent potential given by 

4e 2 

V(r) = f exp(— Pqt 2 ) H erf(/3r), (19) 

where r is the a — a relative distance vector. The parameters of the BFW potential are 
v = -122.6225 MeV, p = 0.22 far 2 , and 8 = 0.75 fin -1 . We choose ^- = 10.4465 
MeVfm 2 and e 2 = 1.44 MeVfm. The three bound states of the BFW potential, ip n em, are 
expanded in terms of the harmonic-oscillator functions V^ m ( r ) = Rn'i{f)Ye m (r) (n' = {q'—£)/2): 



<Pntm(r) = K?R n 'i{r)Y im (r). (20) 

n' 

The phase convention of R n >e(r) is such that it is positive for r greater than the outermost 
nodal point. Table lists the expansion coefficients [/"/. It is seen that the bound state 



TABLE I: Expansion of the three bound states of the BFW potential in terms of the harmonic- 
oscillator functions. See Eq. (|2U|). 
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FIG. 1: Plot of Log(l/7) for N = 1 — 174 eigenvalues obtained for Q max = 30, where 7 is the 
eigenvalue of the operator T defined by Eq. using the three forbidden states of the BFW 
potential (see Eq. 0). In the inset, the eigenvalues around N = 130 are compared with those 
obtained for r H °. 



wave functions (pnem(f) are well approximated by the oscillator functions ipf m (r) (q—2n+£), 
but we will see that the small difference between the two functions produces a huge effect 
on the energy of the 3a system. 

The expansion (J2TJJ) makes it possible to solve accurately the eigenvalue problem (JJJ) in 
the translationally invariant harmonic-oscillator functions provided the maximum number 
of total oscillator quanta Q max is sufficiently large. We treat only + case in this paper. 
For Qmax = 30, there are 444 translationally invariant basis functions, and from them we 
can construct 174 independent basis functions of [3] symmetry. Figure [T] displays Log(l/7) 
, where 7 is the eigenvalue of T obtained for this case. It is clear in the figure that the 
eigenvalues 7 diagonalized in these [3]-symmetry functions are classified into three groups: 
the first group (Group I) is characterized by 7 = and it contains 129 eigenfunctions 
(the values of 7 are actually less than 10~ 14 , which we considered zero within the present 
numerical accuracy). The elements of this group are allowed states and constitute the 
basis for diagonalizing the Hamiltonian in Eq. JTOJ). The second group (Group II) contains 
just two states whose eigenvalues are approximately zero; 71 = 1.35152 x 10~ 5 and 72 = 
1.07152 x 10 -3 . The appearance of these states was suggested in a Faddeev treatment of the 
3a system 0, H^. The third group (Group III) contains 43 elements with eigenvalues 
of order unity (~ 1). These last states must be discarded because they have a substantial 
overlap with the forbidden states. Differently, using r HO , in which the two-body forbidden 
states are the harmonic-oscillator functions as defined in Eq. (|16|). the eigenvalues are divided 
into two groups only corresponding to I and III, and there is no such a group like II. The 
inset of the figure compares the eigenvalues of both cases for N = 126 — 136. 

Table |H] shows the Q max dependence of the number of eigenfunctions belonging to each 
group. There are two noteworthy points; first, the two solutions of Group II always appear 
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TABLE II: The number of allowed and forbidden states for the 3a system with + . The a — a 
forbidden states are taken as the three bound states of the BFW potential. Q max is the maximum 
number of total oscillator quanta, Abasis is the basis dimension of the translationally invariant 
harmonic-oscillator functions with Q = — Qm&x, is the number of totally symmetric functions 
with (5 = — Qmax, -/V 7=0 is the number of allowed states of Group I, A 7 ~° is the number of 
eigenfunctions of Group II, and AH""" 1 is the number of forbidden states of Group III. 
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TABLE III: Expansion coefficients of the two states of Group II in terms of the Pauli-allowed states 
^j3]Q(Afi)«7— rpj ie j nc j ex K j g SU pp r essed when it is unnecessary. See Appendix for the definition 
of the two independent states which appear for Q(Xfi) = 14(64). 
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0.20863 


-0.02453 


14 
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-0.42113 


0.11500 



for all Q m ax values. Moreover, we have confirmed that these solutions (the eigenfunctions 
and the eigenvalues) are very stable against the increase of Q max , e.g., for Qmax = 50, 71 = 
1.35191 x 10~ 5 and 72 = 1.07153 x 10~ 3 . Second, the sum of the members of Groups I and 
II is equal to the number of Pauli-allowed states which occur for the case of = S n r n , 
namely the bound states (f n e m are replaced with the harmonic-oscillator functions ■ 

For later discussion we expand the two states of Group II in terms of the harmonic- 
oscillator basis functions. Table IHII lists the overlap of these two states with some of the 
Pauli-allowed states •^p( A w re7 ~° ( see Appendix for the expressions of these Pauli-allowed 
states). From the table we can notice that the two states with the eigenvalues 71 and 72 
contain predominantly Q = 10, 14 and Q — 8, 12 components, respectively. 
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In order to understand the reason why the orthogonalizing pseudopotential approach 
gives a result quite different from that of the microscopic model, we first diagonalize the 
Hamiltonian H in the following three different basis sets and compare their energies and 
properties for the ground and excited + states. The three sets are (i) a set spanned by 
the basis functions of Group I, (ii) a set spanned by the basis functions of Groups I and 
II, and (iii) a set spanned by the Pauli-allowed states -^-j ) 3 ( j Q ( A ^ K7 ~ . The calculation has 
been performed for some Qmax values to check the convergence of the solution. The results 
obtained with the three sets are listed in Table HV1 The energies of set (i) correspond to the 
eigenvalues of Eq. ()10|). and will approach the energies -E'(A) for A sufficiently large as will be 
shown later. The ground state energy for Qmax = 60 is —0.22 MeV, which is consistent with 
the value obtained by the projection technique [3, The calculated energy is far from 
the experimental value —7.27 MeV and rather close to the 3a threshold where the second + 
state is observed. The root mean square (rms) radius i? rms of the a particle distribution for 
the ground state turns out to be 2.31 fm, which is much larger than 1.92 fm jjjj. The latter 
value is obtained by the 3a boson wave function which is mapped from the microscopic 3a 
wave function that reproduces both the energy and charge radius of the ground state of 12 C. 
The result of set (ii) is in a sharp contrast to case (i). The ground state is now very strongly 
bound and consequently its rms radius is too small compared to the result of the microscopic 
calculation, whereas the calculated second + state appears near the experimental energy 
and its rms radius is increased to about 2.34 fm, though it is still too small compared to 4.26 
fm of the microscopic calculation |l7[ . It is remarkable that the results calculated with sets 
(ii) and (iii) are very similar. From this we may conclude that the function space spanned 
by the members of Groups I and II is nearly the same as that of the Pauli-allowed states. 

As demonstrated above, the two states of Group II play a key role in determining the 
characteristics of the macroscopic 3a model. As shown in Table Hm they contain a significant 
amount of the Pauli-allowed states with Q(Xfx) = 8(04), 10(62), 12(12,0), 14(14,0) etc., 
which should be included from the microscopic point of view. The Hamiltonian matrix 
elements between these two states are given, in units of MeV, as follows: 



We see that the energy of the state with the eigenvalue 72 is already lower than the ground 
state energy obtained in the calculation with set (i). This is due to the fact that it is 
dominated by the lowest shell-model configuration of Q{\jj) =8(04). The energy obtained 
by coupling the two states of Group II is —9.8 MeV, so they clearly contribute to lowering 
the energy. Whether these two states are included in the calculation or not, particularly the 
72 state, is a crucial factor for producing quite different results. 

Next we discuss the origin of the enigmatic behavior of the energy E(\) which is obtained 
in the projection technique. We plot in fig.|2]the energies of the ground and excited + states 
as a function of the projection constant A in the case of Q max = 50. Other choice of 
gives a similar result. The dependence on A follows the pattern observed in Refs. (13. Il5l flif 
which use the basis functions different from the harmonic-oscillator basis. No indication of 
the energy convergence is attained for A < 10 5 MeV, but it begins to show a convergence 
for larger A values. The origin of this behavior is attributed to the particular selection 
procedure of the states in the projection technique. It is clear that the states of Group III 
play no active role from the beginning because the eigenvalue 7 is of order unity and when 
it is multiplied by a large A value the energy expectation value becomes very high. On 




(21) 
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TABLE IV: Comparison of the 
ground and excited + states cal- 
culated in the different basis sets: 
(i) the functions of Group I, (ii) the 
functions of Groups I and II, and 
(hi) the Pauli-allowed states ob- 
tained in the harmonic-oscillator 
basis. E is the energy with respect 
to the 3a threshold. (T), (Vjy), 
and (Vc) denote the expectation 
values of the kinetic energy, the 
nuclear potential energy, and the 
Coulomb potential energy, respec- 
tively. The R vms value of Ref. 0] 
is based on the a boson wave func- 
tion which is obtained from a map- 
ping of the microscopic 3a wave 
function. All energies are in MeV 
and the rms radius is in fm. 
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FIG. 2: Convergence of the energy with respect to the 3a threshold as a function of the projection 
constant A. Q max = 50. 

the contrary, the states of Group I are always active in the diagonalization of the pseudo 
Hamiltonian H as their eigenvalues are zero. The two states of Group II, especially the state 
with the eigenvalue 72 = 1.07153 x 10~ 3 , are in a subtle situation. The expectation value of 
H for this state is 

-7.78 + A 72 , (22) 

so the state is expected to contribute to the ground state significantly as long as A72 is 
smaller than, say 8 MeV, that is, for A < 10 4 MeV, but it will be excluded in the end for 
sufficiently large A. This is confirmed in Table W\ which displays the probability of finding 
the state with 72 in the wave functions of the three + states. The probability found in the 
ground state is 0.64 for A = 10 4 MeV and it decreases rapidly for A > 10 5 MeV. 



TABLE V: A dependence of the probability of finding the state with the eigenvalue 72 in the ground 
and excited + states in 12 C. The calculations are performed for Q max = 50. 
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V. a-a POTENTIAL 



In the previous section we have argued that the two-body forbidden states for the a — a 
relative motion have to be defined referring to the microscopic wave function of the a particle, 
not from the bound states of the a— a potential. Based on this viewpoint, we have shown that 
the BFW potential is so attractive that the ground state energy of 12 C is strongly bound. 
Here we discuss an a — a potential which is more suitable for the macroscopic calculation. 
In doing so, we impose the condition that the a— a relative motion should be orthogonal to 
the 2a Pauli-forbidden states, that is, the harmonic-oscillator functions i/;f m with q = 0,2. 

First we modify the strength of the BFW potential so as to reproduce the + resonance 
energy of the 2a system under the above orthogonality constraint. The resulting potential, 
which we call MBFW, has a strength parameter of vq = —121.888 MeV. Note, however, 
that this modified potential does not reproduce the 2 + resonance energy. Table IVT1 lists the 
energies and -R rms values calculated with this potential in the basis functions of set (iii) with 
Qmax = 50. They are almost the same as those obtained with the BFW potential. Compare 
these with the result of set (iii) in Table HV1 Since no improvement is attained by the MBFW 
potential, we conclude that the change of Vq of the BFW potential does not lead to an a— a 
potential suitable for the macroscopic calculation. 

It seems that the range of the BFW potential has no direct connection to the underlying 
nucleon-nucleon potential or to the wave function of the a particle. Another potential we 
consider here is similar to the a— a folding potential used in Ref. j23[. We assume the same 
form as the BFW potential, and choose the potential parameters as follows: 

v = -107.9 MeV, p = 0.20 far 2 , (3 = 0.589 fin -1 . (23) 

Again, this potential reproduces the + resonance of 8 Be, but predicts the 2 + resonance to 
be about 1.9 MeV, which appears too low compared to the experimental value, E = 3.13 
MeV. The range of this potential is longer than that of the BFW potential, and consequently 
its strength is made considerably weaker. The energies and R rras values of the ground and 
excited + states are obtained with this new potential in the Pauli-allowed space of set 
(iii) and they are listed in Table IVII The energy of the ground state is now pretty much 
improved, though it is still too strongly bound by about 4 MeV. 

The truncation with Q max = 50 may not be good enough to obtain convergent results 
for the excited states. To check this, we diagonalize the Hamiltonian in a-boson functions 
g which are obtained from correlated Gaussians G through the following fermion— >boson 
mapping: 

) = (<Kai)0(a 2 )0(a 3 ) I A{<f>{a x )<j>{a 2 )<l>{a^G{a, R)}), (24) 

with 

G(a, R) = expj - a 12 (R 1 - R 2 f - a 13 {R 1 - R 3 ) 2 - a 23 (R 2 - i* 3 ) 2 }, (25) 

where <fi{a) is the intrinsic wave function of the a particle, and A is an operator which 
antisymmetrizes the 3a microscopic wave function built on the 12-nucleon coordinates. A 
parameter a = (012,013,023) specifies the correlated Gaussian. By choosing the parameter 
a appropriately, one can express a variety of different shapes of the three-body system. 
The function g obviously has [3] symmetry, i.e., g(a,r x , p x ) = g(a,r 2 , p 2 ) = g(a,r 3 , p 3 ) and 
g(a, — 7*1, Pi) = g(a, r x , p 2 ) etc. Furthermore, it is orthogonal to the Pauli-forbidden states 
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TABLE VI: Comparison of the ground and excited + states calculated with the different a — a 
potentials. The basis states are either harmonic-oscillator functions of set (hi) with Q max = 50, or 
the a boson wave functions defined in Eq. (|24|). See the caption of Table ITVl 
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Eq. (231) 
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4.15 


3.17 



i\) q lm (q = 0, 2) if 0(a) is constructed from the (0s) 4 harmonic-oscillator function (with its 
center of mass part being eliminated). Thus the function g, for arbitrary a, satisfies all the 
conditions necessary for the allowed states as discussed in sect. 11111 Moreover, it should be 
noted that the boson function g contains not only the two-body Pauli effects but also full 
three-body exchange effects. The method of calculating g is given in Ref. (l7| . 

The 3a wave function \1/ is now approximated with a combination of g. The parameters 
of g are selected by the stochastic variational method pj. The results of this calculation 
are also given in Table IVTl It is seen that the harmonic-oscillator expansion converges very 
slowly except for the ground state. The R Tms value for the state increases to 3.17 fm, 
which is still smaller than the value calculated in the microscopic model. As noted above, 
the calculation using g includes not only the two-body Pauli effects but also the three-body 
exchange effects, whereas the calculation of set (iii) includes only the two-body Pauli effects. 
The fact that they give virtually the same result for the ground state indicates that the 
three-body exchange effects are negligibly small. Thus we do not agree with the discussion of 
Ref. [lH ,which was made for a possible reason of the sensitivity of the projection technique. 

Though a considerable improvement is obtained with the new potential, the ground state 
energy is still too low by about 4 MeV. One possible reason for this is that the potential 
is too attractive for the d wave (and probably for the g wave as well). Then the lowest 
Pauli-allowed state of Q(Xfi) = 8(04) with L = 0, for example, gains too much energy, 
because the dominating partial waves of this state are d and g waves, as understood from 
the SU(3) coefficients ((40)£, (40)£ || (04)0) 2 , which are 0.28, 0.36, and 0.36 for £ = 0,2, and 
4, respectively. 

The above consideration suggests that a more suitable a — a potential should have an l- 
dependence. In this respect, it is interesting to recall a series of recent publications |26l liil 
3(|, where the utility of a 2a resonating group method (RGM) kernel has been investigated 



for the 3a system. The RGM kernel is non-local and energy- dependent, and therefore a 
solution of the system with such a potential requires some sort of self consistency procedure. 
The ground state energy of 12 C turns out to be quite acceptable, that is, it is not strongly 
bound and is predicted to be around the value obtained by the microscopic calculation. Since 
a non-local potential could be converted to an equivalent local potential with i dependence, 
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an a— a potential suitable for a macroscopic calculation may be found in ^-dependent form. 
VI. SUMMARY 

We studied the ground and excited + states of 12 C in the 3a model, using the deep BFW 
potential. The three redundant states with Os, Is, and Od of this potential were eliminated 
from the function space of the 3a system. Two methods of elimination were tested. One is 
to diagonalize T in a complete set of translationally invariant harmonic-oscillator functions 
with [3] symmetry in order to separate the allowed states from the forbidden ones. The other 
is to use the orthogonalizing pseudopotential Ar and to obtain the energy which converges 
for A large enough to kill the contribution of the forbidden states. These methods, though 
equivalent, were tested in order to clarify not only the origin of the enigmatic A-behavior of 
the second method, but also to give a reason why the ground state energy is very high. 

We found two eigenstates of T with almost zero eigenvalues. One of them especially plays 
a key role. This one has a dominant configuration corresponding to the lowest shell-model 
state of Q(Ayu) = 8(04), and in spite of its small eigenvalue, of order 10~ 3 , it was eliminated as 
a forbidden state in the orthogonalizing pseudopotential approach. The energy expectation 
value of this state alone is already close to the ground state of 12 C. With the exclusion of this 
state, there is no hope to get a ground state energy close to the experiment. In addition, 
we showed that the peculiar dependence of the energy on A can be understood from the 
admixture of this state in the solution: for A < 10 4 MeV, this state can have the effect of 
lowering the energy, but for A > 10 5 MeV it can no longer keep this effect in the ground 
state energy. 

The existence of the states with almost zero eigenvalue of T is not a characteristic of the 
BFW potential, but also applies to other a — a potentials provided their bound states are 
well approximated with the three harmonic-oscillator functions ip q lm (q = 0, 2). We find that 
the potential of Eq. (J25j) produces two almost Pauli-allowed states with ji = 1.59640 x 10~ 5 
and 7 2 = 1.17812 x 10- 3 for g max = 50. 

In our opinion the eigenstates with almost zero eigenvalue should be included in the 
3a allowed space as they correspond to the most important shell-model configurations, 
otherwise, the 3a macroscopic model loses a link to the microscopic 3a model. To include 
these states automatically, it is safe to define the operator F in terms of the harmonic- 
oscillator functions. With this definition of T we demonstrated that the BFW potential 
is too attractive and an appropriate a — a potential has to be determined referring to the 
underlying microscopic model. 
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Appendix: Pauli-allowed states in SU(3) coupled basis 

We give the Pauli-allowed states, Xoo^ 7_ °j for Q = 8 — 14 in terms of the SU(3) 
coupled basis. The index k is suppressed when it is unnecessary. 
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xrr^M^Vi^vjff- (26) 

Q = 10 



+ ^ ! [^ 8 (r 1 )^(p 1 )]i 6 o 2) . (27) 



• g = 12 



x g, 12( oa)^o = ^ (ri) ^ (pi)ff , 

5V5 [*WVi)C + -^[*" , (r I )*V 1 )]S?>. 



4V43 2V43 

Xoo 



gl2(12 ,o h =o = ^ W 4 (ri)v , 8(pi)]S)2 ,o) _ £vl5 [v ,6 (ri ^, (ft)]r , 



116 Lr v /r v ^ i/JUU 58 
+ ^[* 8 (r 1 )»'Aft)]?o 2 ' 0) - g^VO^ftC'"* 

+ ^W 12 (r 1 )^ (p,)]r ) . (28) 

We have two Pauli- allowed states for Q(Xfi) = 14(64), so there is a freedom to fix them. 

Here they are defined in such a way that one of them (« = 2) contains no [V' 10 (^i)V ;4 (Pi)]oci 4 ' ) 
component. 
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xT (26h =° = ^ 6 (n)^V)]of + ^VO^Pi)]^, 
x [3]i4( 6 4) K =i 7 =o = .479702[^(r 1 )^ 10 (p 1 )] 6 4) + 0.619292[^(ri)^ 8 (Pi)C 

+ 0.206431[^(r 1 )^ 6 (Pi)]S 4) +0.586302[^ 10 (r 1 )^ 4 (p 1 )] o 4) , 
x [3]i4(64) K=27 =o = o.337100^ 4 (r 1 )^ 10 (p 1 )] 6 4) - 0.522233[V> 6 (ri)^(Pi)]£ 4) 

+ 0.783349[^ 8 (r 1 )^ 6 (p 1 )] 6 o 4) , 

xS 14(10 ' 2) ^° = p^ir(r^\ Pl )]&°> 2) + M^ 6 (^ 8 (pi)]£°' 2) 



8V19 2V3^ 



8V19 

9\/3 r , S/ \ , fi/ si (14,0) . 21 r/ iQ/ s ,4/ \l(14,0) 



2 v / 742 2V/742 

^wwr + ^[^(ri)^(/h)ir ) . 

(29) 
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